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A new method is proposed for the calculation of full density matrix and thermody- 
namic functions of a many-boson system. Explicit expressions are obtained in the pair 
correlations approximation for an arbitrary temperature. The theory is self-consistent 
in the sence that the calculated properties at low temperatures coincide with that 
of Bogoliubov theory and in the high-temperature limit lead to the results for clas- 
sical non-ideal gas in the random phase approximation. The phase transition is also 
revealed as a concequence of Bose-Einstein condensation deformed by interatomic 
interactions. All the final formulae are written solely via the liquid structure factor 
taken as a source information instead of the interatomic potential and, therefore, in- 
terconnect only observable quantities. This gives also a possibility to study such a 
strongly non-ideal system as liquid 4 He. 
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INTRODUCTION 



The properties of liquid helium are known to have brought about ample literature. Yet due 
to its unique characteristics this quantum fluid keeps attracting the attention of specialists 
in the field of theoretical and experimental physics alike. The first principle microscopic 
description of thermodynamic and structure functions as well as the phenomenon of Bose- 
Einstein condensation at a considerable distance from the absolute zero, in particular in the 
vicinity of the A-transition, is still a problem which cannot be taken as solved to the very 
end. The latter statement can be supported by a fairly instructive study of heat capacity 
in the vicinity of the A-transition. The A-like form of the heat capacity in the vicinity of 
the phase transition of liquid helium into the superfluid state has been taken for a loga- 
rithmic ddivergence with the critical exponent a — > 0. This view found its way into both 
text-books and monographs [H El El IH El EI ■ Precise experiments helped to find out that 
in fact there is no divergence in the heat capacity [H] even though the exponent a is 
indeed a small but negative number, a = —0.01056 [HJ. The studies of the A-transition 
based on the renormalization group method p] provide us with a possibility to carry out a 
correct calculus for the so-called universal characteristics solely, i. e., the critical character- 
istics of the thermodynamic functions and the relations of the amplitudes of their leading 
asymptotics at the temperature on either side tending to the phase transition point. Even 
though the thermodynamic potential functional from the two-component order parameter 
for liquid 4 He was calculated precisely owing to the coherent states depiction OUT], yet 
its subsequent simplifications necessary for the implementation of the renormalized group 
approach make it impossible to describe the system's characteristics outside the closest 
vicinity of the phase transition point using the same method. Notwithstanding tangible 
efforts of the researchers the renormalization group method did not yield the logarithmic 
divergence of the heat capacity (the a exponent was received as a small but still finite pos- 
itive number, the power divergence having been obtained). Only in the subsequent studies 
which made use of the summation procedure of the Borel perturbation theory divergent 
series established the negative value of the exponent: a = —0.01294 [12], a = —0.0150 [13J, 
a = -0.01126 [H]. 

The fact that heat capacity at the A-point acquires a finite value finds itself in agreement 
with F. London [TH] in the sense that the A-transition in liquid 4 He is a Bose-condensation 
deformed by the interatomic interaction inherent of the ideal Bose-gas bringing about the 
steepness of the heat capacity curve in the vicinity of the Bose-condensation point. It is not 
easy to reveal the connection between superfluidity and Bose-condensation unequivocally. 
This connection is a lot more complex than a simple correspondence. For instance, at the 
temperature equaling absolute zero all the atoms of the ideal Bose-gas find themselves in 
the states with a zero momentum thus forming the one-hundred per cent Bose-condensate 
(BC). However, such a system is not superfluid. Conversely, liquid 4 He at the absolute 
temperature zero is superfluid even though the number of atoms in the BC, as showed 
both by theoretical and experimental studies, is but a small part of the entire number of 
atoms. It is also well-known that one- and two-dimensional systems tend to reveal superfluid 
characteristics with the BC being absent. Here we are concerned with the rise of the so- 
called non-diagonal long-range order when the one-particle density matrix in the coordinates 
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description dwindles in compliance with the power law rather than the exponent law, at the 
matrix arguments being disentangled over infinity. On the other hand a mere assumption 
that the BC exists in the weakly imperfect Bose-gas made it possible for Bogoliubov [16j 
to obtain from the first principles the energy spectrum whose properties are close to those 
of the liquid 4 He. 

One other fairly complex task lies in according Bogoliubov's theory [[16j well-suited for 
the low temperatures with the random phase approximation in the theory of classical sys- 
tems. It may be possible to find expressions for the thermodynamic and structure functions 
of the Bose-liquid which would yield the results of Bogoliubov's theory for the temperatures 
T^OK. For high temperatures in the quasi-classical h — > limit this expression would 
bring closer to the random phase approximation of the classical non-ideal gas theory. It 
might be supposed that such an expression would also give good results in the intermedi- 
ate temperature area where the transition point of liquid 4 He into the superfluid state is 
located. It is clear that we are referring not to the general formulae at the level of defi- 
nitions which just cannot "be brought to a certain number", but to the discovery of such 
a first principles method of calculating the thermodynamic properties of the many-boson 
system which starts from the N particles Hamiltonian and makes it possible to put forward 
a regular perturbation theory with a consistent consideration of many-particle correlations 
starting with the two-particle correlations. 

In the present paper we mean to dwell upon one of the possibilities of solving these 
problems. Some of the results to be presented here were briefly summarized in [T2| just as 
an illustration of certain tricks applicable to the wave functions and the statistical operator 
[18j. However, the full scope of these findings has never been published. We shall calculate 
explicitly the complete density matrix taking into account the two-particletwo-particle 
interatomic correlations as well as the Helmholz free energy, the energy and also a pair 
liquid structure factor which will abide by the above requirements. 

1 INITIAL EQUATIONS 

Let us consider the N set of spinless Bose-particles with the mass m with the Hamiltonian 

N ~2 
j=l l<i<j<N 

where the first term is the operator of the kinetic energy, pj is the operator of the momen- 
tum of the j-th particle. The second term presents the potential energy of the two-particle 
interaction between the particles with the coordinates x = (ri, rjv). The particles move- 
ment is limited within a volume of V. 

Let ip n (x) be taken for a system of eigenfunctions of the Hamiltonian H and E n are its 
eigenvalues. Let us consider such an equation for the statistical operator: 

e-^M*) = e~ PEn Mx), (1-2) 
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where (3 = 1/T, T is the temperature of the considered system of particles. Let us introduce 
a certain arbitrary function if = <p(x) and let us multiply the left-hand side of equation 
(1.2) by (p(x). We will integrate it over all coordinates x: 

J ^{x)e~ p " i) n {x)dx = e~ pEn J ip(x)ip n (x)dx, (1.3) 

where 

J dx = J dr 1 ...J dr N . 
We will impose the following condition for the function (p(x): 

<p(x)4> n (x)dx ± 0. (1.4) 



Then making use of the self-conjugation of the operator H we shall transfer it in equation 
(1.3) from the function i/) n (x) to <f(x) multiplying the 1. h. s. side of this equation by ipn(x') 
and summing it by all the indices of the n states: 

]T / r n (x')il>n(x)e- f>6 <p(x)dx= f (pW^nWe-^M^dx. 

n n 

We make use of the completeness condition 

i>*n{x')i>n{x) = 5(x' - X), 

n 

and also proceeding from the previous equation we will obtain the following equation with 
the permutation of the variables x for x' and inversely: 

e~ p6 (p{x) = J (fi(x')R N (x'\x)dx f , (1.5) 

where 

R N (x'\x) = ^C(^)e- /3£ ^n(x') (1.6) 

n 

is the density matrix in the coordinate representation. 

In true fact equation (1.5) is self-evident and could have been written out at once as a 
matrix presentation of the statistical operator action over any function (p(x). In the Dirac 
notations this equation can be written as follows 

(x\e- p6 (p) = J (x\e-^\x')(x'\ip)dx\ 

(1.7) 

(x\e-P"\x') = R N (x'\x), (x'\<p) = (j)(x'). 
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It is clear that formula (1.5) or (1.7) is applicable both to one-particle and two-particles 
problems. Notwithstanding the fact that it is obvious, this formula makes it possible to 
fully restore the density matrix (1.6). Here most significant heuristically is an arbitrary 
choice of the function (p(x) on condition that it is non-orthogonal to the eigen-functions of 
the H operator and the existence of the integral (1.4). 

The present article is concerned with the calculus of the Rn(x'\x) density matrix for 
the studied many-boson system from the initial equation (1.5). 

Should the interparticle interaction be switched off, the density matrix will transfer 
into that of the ideal Bose-particles N system. That is why we will try to construct the 
density matrix of the N system of interacting particles as a product of the density matrix 
of the ideal Bose-gas R° N (x'\x) and the P^(x'\x) factor taking into account the interparticle 
interaction: 

R N (x'\x) = R° N (x'\x)P N (x'\x), (1.8) 

where at temperature T^O 

3N/2 

N\ \27if3h 2 



I f m* \ 



* N 

m 



2(3h 2 

j 



(1.9) 



the summation over Q is at the same time the summation over all the permutation iV! 
numbering the particles coordinates; at T = 0, when all the particles have zero momenta, 
the matrix R%(x\x') = 1/V N . The mass of the particle in (1.9) is understood as a certain 
effective mass m* which equals the initial mass m unless the interaction takes place. Thus 
a part of the contribution from the interparticle interactions is taken into account by 
renormalizing the particle mass, the remainder being left in the P^(x'\x) factor. Postulating 
this for the density matrix (1.8) is justified also by the fact that it was obtained under certain 
assumptions by a different method in jTH] by a direct calculus from (1.6). The expression 
for the effective mass was found to be in agreement with the formula for the effective mass 
of the impurity atom moving in the Bose-liquid if the initial impurity mass coincides with 
the mass of the atom of liquid. Yet for the R^{x'\x) matrix in [THj an expression was found 
which is correct only at T — > K. Here we will try to find a formula for P^(x'\x) which 
will be capable of working in the entire temperature interval. We will also postpone the 
discussion of the issue of determining the effective mass contained in expression (1.9). 

Falling back on the phenomenological consideration concerning the insignificance of the 
atoms permutation which are located at the distances smaller than their own sizes quite 
some time ago R. Feynman (2H1 E| constructed the iV-particle distribution function of the 
liquid 4 He, i. e., the diagonal elements of the density matrix as a product of the ideal Bose- 
gas distribution function multiplied by the factor accounting for the atom impenetrability. 
We would like to base our analysis on the precise equation (1.5) and suggest a consistent 
method for the calculus of the Pj^{x'\x) function computing the latter explicitly accounting 
for the two-particle correlations. We will choose the Pn(x'\x) matrix which takes into 
account the interparticle correlations as follows: 

P N {x'\x) = exp jc + J^ci(g)pqp_ q - ^^c 2 (g) p q P- q + p q p'_ q \, (1.10) 

^ q^O q^O J 



5 



where the Fourier coefficients of the particles density fluctuations 

1 N 

iV 



at q 7^ 0. The components of the wave vector q cover the integer values devisible by 
2n/V 1 / D where D is dimentionality of the cubic box which contains the system of particles 
studied. We will determine the coefficient functions Co, ci(q), 02(g) from equation (1.5). 

We have confined ourselves to the consideration of two-particle interparticle correlations 
in expression (1.10). Accounting of the three-particle and higher correlations is accomplish- 
able by adding the members with a product of three or more p q in the exponent in (1.10). 
In our work we will not take them into account explicitly, yet we will return to them when 
discussing the issue of the effective mass m* . The exponent form of the Pn(x'\x) matrix is 
caused by the classical boundary of its diagonal elements when Pn(x\x) turns into the usual 
Boltzman factor exp(— /?<&), where <3> is the potential energy of the interparticle interaction 
thus equaling the second term in the Hamiltonian (1.1). 

Finally, we must also choose the appropriate ip function which is contained in equation 
(1.5) and meets requirement (1.4). In our case we will adjust it from the class of such 
functions: 



(p(x) = exp 



E A (^ 

Lq^O 



(1.13) 



where the arbitrary coefficient function X(q) is a real function from the wave vector module 
q — |q| with all the necessary properties. 



2 THE LEFT-HAND SIDE OF EQUATION (1.5) 

Let us represent the action of the statistical operator on the p function in (1.5) in the form 
of the exponent function: 

e u = e^V (2.1) 

Here again when fixing the exponent form of this function we proceed from the assumption 
that in the quasi-classical limit when K — > the statistical operator comes up to the 
Boltzman factor which just as the (p function has an exponent form. Differentiating either 
sides of the equation with respect to inverse temperature (3 we will find the equation for 
the unknown function U : 

dp 
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or explicitly 



_ 2 N 



~% = -^£[^+W> 2 ]+*. (2-2) 
where from (2.1) and (1.13) it is apparent that 

U = Y,K<l)P<v when £ = 0- (2-3) 

q^O 

To solve equation (2.2) we will choose the function U in the form of the series "by the 
degrees"of the p q quantities: 

U = a + ^ a l(?)Pq + \Yl a 2(<?)PqP-q + • • • ■ (2.4) 
q^O q^O 

We have agreed to take into account the two-particle interparticle correlation only that is 
why the series is truncated at the second degree of p q . Higher correlations denoted as dots 
in (2.4) can be taken into account by summing higher degrees of p q . It is convenient to 
work with our equation in terms of the p q quantity rather than the individual coordinates 
x. That is why we will show the potential energy <3> through the quantities of p q : 

N(N - 1) N , . , . 

* = \y V* + ^7 E ^(A'qP-q " X )> ( 2 ^) 

q^O 

where the Fourier coefficients of the two particles potential energy 

v q = J e- tqR $(i?)rfR. (2.6) 
Equation (2.2) will now look as follows 

du _ v h2k2 ( du ° 2u du du 

d[3 2m Y «9pk <9p k <9p_ k <9p k <9p_ k 

1 ^ ^ fi(kk') _ ( d 2 U t dU dU 

k+kVO 

N(N-l) ^ N , . , , 

+ 2y + 2^ 2y^(PkP-k - 1). (2.7) 

k^O 

We substitute in this equation expression (2.4) for the U function and from the equivalence 
condition of the coefficient functions at the equal degrees of p k in the left- and right-hand 
sides of this equation and again taking into account just the two-particle correlation we 
find a system of three equations for the unknown quantities do, and 02(g) 

da v^Vr 2 , N(N-l) ^ N , . 

^ q^O q^O 
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da 'M- fiV - 0l (,)[l-2a 2 ( 9) ], (2.9) 



dp 2m 

dci2(q) h 2 q 2 r . . 9 . N . . 

■Kg) -a 2 2 (q)\ + -u q . (2.10) 



d/3 m ZV ^ J 1/ 

We begin to solve the system of these equations from equation (2.10) for a,2(q). Notwith- 
standing the fact that this equation is nonlinear its solution can be found analytically which 
was suggested for the first time in [21 J in the framework of the so-called shifts and collective 
variables method: 

oc q - 1 1 - e" 2 ^) , . 

^ = - 2 i + a zl e -**(,) ' ^ 




where 



a q = ^|l + —v n I (2.12) 



h 2 a 2 

E{q)=aq ^m~ (2 " 13) 



is Bogoliubov's elementary excitation spectrum [TH] . 

After this we solve in turn equations (2.9) and (2.8). Fairly simple even though cumber- 
some transformations taking into account the "initial"condition (2.3) lead to the following 
result 



2a a 



-PE{q) 



aAq) = X(q) " g+1 1 , (2.14) 



fl 



I £ In 



a q +l 

Oin—1 



2 



1 + a^ 1 c -20E(a) 
a q +l 

1 



where 



+ E A (?) r — ^1 77p (2-15) 

^ (a, + 1) [l + Sie-Wri 

N(N-l) v-^-V, x2 

q^O 

is the energy of the ground state of the many-boson system in Bogoliubov's approximation 

HB|. 

Thus in the approximation of pair correlation we have found explicitly the result of 
the action of statistical operators on the ip function, i.e., the left-hand side of our main 
equation (1.5) which taking into account (2.1) looks as follows: 



e u = J <f(x')R(x'\x)dx'. (2.17) 
To go further we must calculate the integral over x' in the right-hand side of this equation. 
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3 THE RIGHT-HAND SIDE OF EQUATION (1.5) 



Let us now proceed to the finding of the explicit form of the right-hand side of our main 
equation (2.17) making use of the postulated form of the (1.8)— (1.10) statistical operator. 
We have 



ip(x')R(x'\x)dx' 



... R° N (x'\x) expl Y A (<?)Pq + c + Y c i(<?)PqP-q 

J J Wo q^O 



Y PqP-q + ^/-q r dr>1 ■ ■ ■ dr ' N ' 



q^O 



(3.1) 



Let us pass over in this expression from the integration of particles r' l5 by the individ- 
ual coordinates to the integration over the p q variables which are determined by formula 
(1.12). Such a transition is performed by means of Zubarev's [22] transition function which 
is a product of the Dirac 5-function and equation (3.1) and looks as follows: 



'cp(x')R(x'\x)dx' = /•••/ exp<^ Y %)p q + c o + Y p '<iP- c 

J J WO q^O 

\Y C2 ^P^ p -^ + P'qP-q} f Up')( d p'), 



q^O 



(3.2) 



where the weight function 

Up'] 



J...J r° n (x'\x) n' 5 U - -j= E e_iqr M d ^ ■ ■ ■ dr N- 

q^O V V j = l J 



An element of the p' space volume 



(3.3) 



(3.4) 



q^O 



where p q , p q are a real and respectively an imaginary part of the variable p q = p q — zp q . 
The integration by p q , p q takes place in the infinite volume. The prime (') at the product 
sign in (3.3) and (3.4) means that the values of the wave vector q are considered only from 
the half-space of all of its possible values because there exists the dependence p q = p_ q ' or 

Pq = P-q, Pq = "P-q- 

We will use the integral representation for the (^-function and write the transition func- 
tion (3.3) as 



Up') 



. . . j R%(x'\x) I (don) exp 
(dui) exp 



mY^Pq 
q^o 



/ 1 N 

q^O \ ViV j = l 

1 N 

-™ E 77X7^ E « 



-?qr 



«qr. 



q^O 



i=i 



dr[ . . . dr' N 

(3.5) 



here the element of the uj space 



where Wq, o; q are a real and respectively an imaginary part of the u; q complex variable linked 
to p q , the angle brackets denoting the integration over the individual primed coordinates 
with the weight equaling the density matrix of the ideal gas: 



(...) 



.fl N (x'\x)dr' 1 ...dr' } 



N- 



(3.6) 



The integral from the R%(x'\x) density matrix itself by x' equals unity. Indeed, 
j ... j R%{x'\x)dv' l ...dv' N 



1 / m* 
M \2n(3h 2 



3N/2 



E 



dr[ . . . dr' N 

1 / m* \ 3N/2 ^ f f ™* V iV R 2 



(3.7) 



where we have passed to the new variables of integration 

R j = r 'j ~ r Qi, 

and at the same time we assume that the size of the box where the system of particles is 
located is already infinite, V — > oo. We will refer to operation (3.6) as to averaging. 

Let us now pass over to the calculation of the weight function (3.5). We will reflect 
the average from the exponent in (3.5) in the form of the exponent from the non-reducible 
means limiting ourselves to two-particle correlations again: 



Up) = f(d^)e^o^ ew f ni J2^(f2^ 

J < q^O V7V \j = l 



+ EE 



(7T^ qi )(7r^q 2 ) 



qi^0q 2 ^0 

N 

jqir 



- (E« 

\ii=i 



n 



N 

N 

E 

\j2 = l 



N 



N 



E^" ir; 'E 

\h=i 32=1 



(3.8) 



+ 



Let us calculate the averege quantities necessary for us. Thus similarly to (3.7) we have 



N 

E' 

j=i 



-«qr. 



1 / m* 



3N/2 



N 



EE 

Q 1=1 



dr[ . . . dr' N e~ i(i <e~^^ l{r 'i~ rQj? 
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* \ 3/V/2 n r r n 

m V 2^J E E ^ TQl J ■JdR 1 ... dR N e-^e-^^ R " 



" AT! \2ir{3h 2 J q{~{\ m * J \ m* J { } 



Incidentally, it is not difficult to show that the mean for the many-fermion system equals 
zero. This is connected with the fact that the R%(x'\x) density matrix for the fermions is 
antisymmetric as regards the permutation of the primed (r^, . . . , r' N ) coordinates, the value 
of (1.12) being symmetrical. Otherwise said, for Fermi statistics the ip function from (1.13) 
does not meet condition (1.4). Thus, it is necessary to find another function which would 
be non-orthogonal to the wave functions of the many-fermion system. The resuls of the 
calculation of the thermodynamic and structure functions for the Fermi particles on the 
basis of the suggested approach will be published separately. 

Our further steps are analogous to those in the case of (3.9): 

N N N N N 

^2 e~ iqir ^ ^2 e ~ imr ' n ) = ( e~ i(qi+q2)r ^i ) + E E ^e" iqir ^~ iq2r ^ 

jl=l ]2=1 3\=1 31=1 32=1 

2 2 N N N Jh 2 <$ ^19 



= e ~P ^ e" i(qi+q2)rjl + e " QirJl E e ~ iq2rJ2e ^ ' ( 3 - 10 ) 

jl = l jl=l J2=l 

2 N 

ii=i 

Now, taking into account (3.9) and (3.10) we will find that 

TV N N N 

e~ iqir - E e ~ imr ' n ) - ( E e " qir;i )( E e ~ m2r ' n ) ( 3 - n ) 

ii=l 32=1 31=1 32=1 

= _ e -^ h ^ NS ( qi + q 2 ) + [ e -^( qi+q2 ) 2 - e~^^ + ^] >/JVp qi+q3 . 

Here, the second term will be omitted in the accepted approximation of two-particle corre- 
lations when substituting (3.11) in (3.8). Expressions (3.9), (3.11) give correct asymptotics 
at absolute zero temperature, (3 — > oo, as well. 

Further, the weight function in (3.8) can be calculated easily as the integrals by o>£ and 
u;£ are reducible to the Poisson integrals: 



Up') 



j (du) exp ~ \™ q \ 2 (l - e W 2 ™* ) + ™ E u « \P'i ~ p * e 

■' q#0 q#0 
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cxp 



2 ^ 

q^O 



B 2 2 |2 



1 — e P 2m* 



II / _ 2 /3^2 2 _\ 
q^O 7T [ 1 - e ^2-»* J 



Now, just as in (3.12) we can carry out integration in (3.2): 



(3.12) 



/ <p(x')R(x'\x)dx' = exp < c H — j 



A 2 (g) 



x 



exp <^ ^A(g)p q 



x exp 



[ft (9) 



2 ^ 



2 ^ 

q^O 

c 2 (g) + 



c i(e) 



1 — e P 2m* 



C 2(<?) + 1 _ e -2^2 9 2 /2m , 

c 2 (g) 



1 - e~ 2 ^ 



^2 2 / ^2 2 

+ e-^/(l-e- 2 ^' 



c 2 (g) + l/(l-e- 2 ^)] 



PqP-q 



e~ 2 « 



1 - e" 2 ' 3 2 



fc 2 q 2 



PqP- 



ir 

q^O 



1 + (1 -e-2^V/2m*) C2 ( g ) 



. (3.13) 



Thus, we have also found the right-hand side of the main equation (1.15). Now we can 
proceed to determining the unknown coefficient functions Co, Ci(q) and 02(g). 



4 EQUATIONS FOR THE c n (q) COEFFICIENT 
FUNCTIONS 

To satisfy condition (2.17) the exponent factor in (3.3) should equal the U function from 
(2.4) for any values of the p q variable and the arbitrary function X(q). From this condition 
we find the equation for the unknown fucnctions cq, ci(q), C2(q) . Thus, we will equate the 
coefficients at the identical degrees of p q in (2.4) as well as in the exponent of the right-hand 
side part of equation (3.13): 



a 



ai(g) 
02(g) 



= co + ^A 2 (,) 



q^O 



c 2 (g) + 



1 



2„2 

1 _ e - 2 ^ 



\ In [l+(l- e - 2 ^ V/2m *)c 2 (g)], 

q^O 



A(?) 



ci(q) + 



^ 2m* 



1 



Z P 2m* 



Cl(<?) + 



1 - e -2/« 



02(g) + 

c 2 (g) + 



1 - e 
1 



-2/34 



1 - e" 2 ^ 



(4.1) 
(4.2) 
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c 2 (g) + 



1 - e^- 



(4.3) 



Let us address equation (4.1) taking into account expression (2.5) for ao as a result of 
the arbitrariness of the function X(q). We will obtain one more additional equation equaling 
the multipliers near A 2 (g) on either side in equation (4.1). In equation (4.2) the value of \(q) 
falls out as can be seen from (2.4). Thus, we now have three equations for two coefficients 
Ci(q) and c 2 (g). If our theory is consistent, one of the equations from system (4.1)-(4.3) 
will be satified identically: 



2a q ( 
a„+l 



-pE{q) 



C l(g) 
P2(9)' 



1 _|_ ^r± e -20E(g) 



a q +l 



C2(q) 



-2f3E(q) 



a q + 1 1 + 2azi e -2/H5(,) 2c 2 (g)' 



a 9 +l 

here we have introduced abbreviated notations 



(4.4) 



ci(g) 
02(9) 



ci(g) 
c 2 (g) 



-2/34 



1 _ e -2/3^ 



From the third equation of system (4.4) we have 



c 2 (g) = 



1 , ?2Z± e -2PE( q ) 

a q + 1 1 + q 9+ i e 

2 1 - e -2/3B(9) 

1 + a g coth[/?£(g)] 



(4.5) 
(4.6) 



and from (4.6) 



c 2(<?) = 2 1 a g coth[ / 9£'(g)] - coth 



P 



k 2 q 2 
2m* 



Now, from the first equation (4.4) we find 
and from (4.5) we have 



2 smh[PE(q)Y 



(4.7) 



sinh[(3E(q)\ ginh 



t J 2m* 



(4.8) 
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We can easily check now that the second equation from system (4.4) is indeed met 
identically as expected. 

Finally, from (4.1) we also find the coefficient 



c 



-PE + 2 E ln I _ e -2^ (g) + 2 E ln a <- 

q^O V / q^O 



(4.9) 



5 DENSITY MATRIX AND THE PARTITION 
FUNCTION 

With the help of expressions (4.7)-(4.9) taking into account (1.8) (1.10) we find the explicit 
expression for the iV-particle density matrix: 



R N (x'\x) = R° N (x'\x) expl-pE + - J^ln 



q^O 



a q tanh [|_E(g)j \ 



1 -e 



x 



tanh 
h 2 q 2 



H 4m* 



,2 2 \ i 
' 2m* \ 



-pE{q) 



ex p(-i J2 ( a i coth \P E b)\ - coth ) 

I q^O ^ L J / 



(PqP-q + PqP-q) 



(5.1) 



1 \- 



q^O 



sinh sinh 



fJ 2m* 



PaP- 



By integrating the density matrix diagonal elements in compliance with (1.6) we obtain 
the partition function 



Zn = J Rn{x\x)(1x, 



(5.2) 



from which we will find free energy F = —T ln Z^, and from the latter other thermodynamic 
functions. 

In (5.2) let us pass over from the integration over the individual coordinates (ri, . . . , rjv) 
to the integration over the collective variables (1.11) similarly to what was done in (3.1): 



1 



Z N = Z° N exp\ -(3E + - J]ln 



a tanh [§£(?)] \ 



+ 5> rr 



„2 2 

1 — e P 2m* 



-(3E(q) 



q^O 



X 



exp < — - (a tanh 



q^O 



<,,u \ tanh 
\dp)J{p) 





l J Am* 



-E{q) 



— tanh 







hrq 
4m* 



PqP-q 



(5.3) 
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where the weight function 



tion / , / i N N \ 

'M = II/^-TfE^) >. 

q^O v ?=1 



(5.4) 



q^O • " ' J = 

here the brackets marked with "0" denote the following averaging: 

1 



<...>' 



.)R° N (x\x)dx, 



and the partition function of the ideal gas 

Z° N = R%(x\x)dx. 



(5.5) 



(5.6) 



We calculate the J(p) weight function from (5.4) in the same way as the Jo(p') function 
from (3.3), clearly just with a different averaging operation (5.5) instead of (3.6). Thus, 
by the respective permutation for J(p) we have expressions (3. 5) and (3.8). The necessary 
averages values contained in (3.8) can now be calculated easily as well 



1 N 



and 



1 N 



jl=l 



1 N 



,-iq 2 r.,-2 



32=1 



5 , o(gi)*(qi + q 2 ), Qi^O, q 2 ^ 0, (5.7) 



where S (q) by definition is a pair structure factor of the ideal Bose-gas. The expression 
for it is well-known 



S (q) = 1 + ^n„n 



|p+q|> 



(5.8) 



where 



n. 



-1 g h ' 2 P 2 -i 

Z Q e P 2m* — 1 



(5.9) 



is an average number of particles whose momentum equals hp, the activity of the ideal gas 
z being excluded from the condition 



(5.10) 



Let us also remember the expression for the partition function of the ideal Bose-gas (5.6) 
the particles of which have the mass m* 



N = eX P 



-^ln(l-^ e-^ 2 / 2m *) 



(5.11) 
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The necessary integral by in (3.8) taking into account pair correlations can be easily 
taken. We will find the wave function (5.4): 



m = n' 



1 



vrS'o(g) 



exp 



4E 



2 ^ S (q) 



(5.12) 



Now, taking into account (5.12) the integration over the variables p q in (5.3) is reduced 
to the calculation of the Poisson integral: 



J N 



Z° N exp\-pE + J2 ln 



zE ln 



q^O 



1 +So(q) ( a q tanh 




tanh [§E(q)] 



tanh 



t9£_ 

Am* 



(5.13) 



Expressions (5.1) and (5.13) for the density matrix and for the statistical operator are 
the starting formulae for the calculation of the thermodynamic and structure functions of 
the Bose-liquid. Before we pass over to the calculation of these values we will make some 
preliminary studies. 

First of all we will say that we have not obtained the equation for the ideal mass m*. 
It is not surprising at all as we were working in the approximation of one sum by the wave 
vector q and the difference between the effective mass m* and the initial particle mass m 
is proportional to the sum by q. 

Let us pass to the discussion of the obtained expressions for the partition function and 
the density matrix. It is apparent that when we switch off the interparticle interaction with 



v q = 0, i. e. 



a g = 1 we will obtain from (5.1) that for any temperature 



Rn(x\x') = R 



N 



and 



7° 
N' 



Notwithstanding the naturalness of this condition it will not be met, for instance, by the 
well-known Penrose formula for the iV-particles density matrix. This formula was also 
obtained by E. Feenberg [22]. Applying the method of coherent states it was also found in 
[24]. With the help of wave functions of the many-boson Bogoliubov-Zubarev system [26J 
it was calculated in [23 EH]- Using our notations it looks as follows: 



Rm{x\x') 



n' a « tanh 4f I exp { ~ y ~~ s in t 1 

-q^O J L q^O 



-E(q)/T 



) 



qP-q + PJ 



qr-qy 



q^O 



4 ^ smh[E(q)/T] 



[(PqP-q + PqP- q ) 



x cosh 



E{q) 
T 



(PqP-q + PqP-q)] |- 



(5.14) 
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This expression holds true only for T — ► 0. 
expression for the ideal gas density matrix 



If a q = 1, then from (5.14) we obtain an 



R%(x\x>) = J]' tanh exp { In (l - e~^^) 



q^O 



+ J 2^(^-1 + P'<*P'-J ~ 4 22 si nh(^g 2 /2mT) 



q^O 



[(PqP-q + PqP-q) COSh 



k 2 q 2 
2mT 



(PqP-q + PqP-c 



(5.15) 



which does not coincide with the precise expression (1.9) for m* = m. It is not surprising 
at all as in the p q representation the kinetic energy operator is not diagonal (which can be 
shown at least from equation (2.7)). It is just on its off-diagonal part that we are building the 
theory of perturbations in which each subsequent member of the series has in comparison 
with the previous member an additional summation by the wave vector q. 

It is curious that formulae (5.14) and (5.15) prompt to us how one can formally obtain 
our result for the density matrix. Thus we have the equation 



R N (x\x') = R° N (x\x r 



-R7v(x|a/)(from formula (5.14)) 
R%(x\x')(from formula (5.15)) 



(5.16) 



If we use expressions (5.14) and (5.15) in the right-hand part of this equation for the 
density matrices relation, we will immediately arrive at formula (5.1) with m* = m. This 
trick also hints at a possibility of calculating the effective mass m* . Thus, should we find 
the following approximation for the density matrix (5.14), for instance by solving Bloch's 
equation directly, as suggested in (2HJ , the result will be the following. Firstly, we will 
obtain in the (5.14) exponent the ~ PqiPq 2 Pq3 members from q3 = — qi — q2, viz., a 
contribution from the three-particle correlations and also a share of the contribution from 
the four-particle correlations ~ PqiP-qiPq 2 P-q 2 (appropriately systematized as regards the 
primed and the nonprimed variables (1.11) and (1.12)). Besides, there also arise corrections 
~ X}q/o(- • •) / 'N to the zero approximation coefficient functions around p q p~ q in (5.14) 
(again systematized as regards the primed and the nonprimed variables). A part of these 
corrections can be "hidden"quite naturally by renormalizing the particle mass m. After this 
fixing of the renormalized mass m* we suggest that a q — > 1 and find the ideal gas density 
matrix R%(x\x') in which instead of m we will have m* . Then we will address formula 
(5.16) and by the same reasoning we arrive at (5.1) with the known value of m* . Now from 
formula (5.16) we can also answer the question why formula (5.1) holds true for any mass 
m* . The matter is that me multiply and divide by the same value R° N (x\x') "finding"the 
denominator for all the orders of the perturbation theory by the sums number and by the 
wave vector. 

In the approach suggested here the effective mass also arises in the natural fashion if we 
take into consideration the contribution of the many particle corrrelations in expressions 
(1.10) and (2.4) and if we renormalize the mass by the contributions to be factorized 



17 



(possibly with the dependence upon the wave vector q)). We will not discuss here higher 
approximations and the discussion of the issue of m* will be returned to in Section VIII. 

If we direct the temperature towards zero, i. e., (3 — ► oo then from (5.1) with the consid- 
eration of the fact that the density matrix of the ideal Bose-gas, which in this case is fully 
degenerate, equals 1/V N , we will obtain in compliance with designation (1.6) the following 
expression 



Rn(x\x') 



where the normalized wave function of the main state of the interacting Bose-particles 
system 



fyW 



/a, 



) exp "iX^ - tfPiP- 



coincides with that discovered for the first time by Bogoliubov and Zubarev [2f 

Finally we will consider the classical limit h — > of the density matrix diagonal elements 
(5.1) when p q = p q . The energy E from (2.16) in this area boundary transforms into 



E 



N(N-l) N 



(5.17) 



Then, the first logarithm as a sum of q in the exponent in (5.1) (together with the factor 
1/2) is reduced at K — > to In a q and the second one to ln(l/a g ); they cancel each other, 
the factor at p q p- q equaling 



a q tanh 



/3 



E(q) 



tanh 







H 2 q 2 
4m 



(5.18) 



We have made use of the fact that at high temperatures m* 
from (1.9) at K — ► only an identical permutation survives 



m. In the R%(x\x) matrix 



1 



m 



N\ \2-k0W 



3N/2 



Consequently, bringing everything together we have the exact classical solution for the 
density matrix 

3N/2 



Rm(x\x) 



(. 



m 



N\ \2ir(3h 2 , 

where the potential energy $ is given by equation (2.5). 
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For the partition function(5.13) taking into account the fact that in the classical limit 
h — > the structural factor of the ideal gas (5.8) So(q) = 1 we find a well-known expression 
in the random phase approximation j3| I25|: 



Z° N exp { - i 



N(N 



2V 



7° 



111 

yN 

77T 



1+Py"* 



-v 



R N 



(5.19) 



m 



( 2tt/3^ 2 



3./V/2 



Thus from expressions (5,1), (5,13) for the density matrix and the partition function 
we have received all the known limiting cases both in the essentially quantum and classical 
regions. For that matter we can expect these expressions to give good results also in the 
intermediate temperature region, in particular in the vicinity of the A-transition point. 



6 ENERGY 



From expression (5.13) for the partition function we find the energy by the well-known 
thermodynamic equation: 

E =ik (,3F) = -ik nZN - 



Simple calculations give 

E : 



+ 



x 



+ 



v ^ E(q) h 2 q 2 

E ° + 2^ p f3E(q) _ 1 + 1^ 



q^O q^O 



i ,3/i 2 9 2 
Z e 2m* — l 



e P 2m* — 1 



- i£ 



q^O 



E{q) 



H 2 q 2 /2m** 



1 E- 

2 1 

q^O 



So(q) 



sinh[f3E(q)\ sinh[/3/?, 2 g 2 /2m*] 

1 



+ S (q) (a^tanh [f3E(q)/2] - tanh [f3h 2 q 2 / Am*]) 



a q E(q) 



H 2 q 2 /2m** 



2 [cosh 2 [f3E(q)/ 2} cosh 2 [(3h 2 q 2 /Am*] . 



dS (q) 

d/3 



a q tanh 



'§*<«) 



tanh 



k 2 q 2 
Am* 



(6.1) 



where the value 



m = m 



1 + (3m* 



d_ 

d(3 \m 



(6.2) 
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For the actual calculus of the temperature dependence of the energy E we will write in 
more detail the ideal Bose-gas structural factor (5.8) and its derivative with respect to the 
inverse temperature 



dS (q) 

d(3 



2 ^ 9n p 
p 



n |p+q|> 



(6.3) 



which enter (6.1) with the consideration of the Bose-Einstein condensation phenomenon. 

Bose-condensation arises on condition that zq = 1 and the critical temperature T c is 
determined from equation (5.10) jTHHIH] 



Tr 



2irh 2 



vnr 



N 
V, 



am 



1 2 /3 



(6.4) 



where Riemann ^-function C(3/2) = 2.612375... . As the effective mass also happens to be a 
function of temperature, expression (6.4) as a matter of fact is the equation for determining 
the temperature of the Bose-condensation T c . 

In (5.4) we will single out the terms with the average number of particles no whose 
momenta equal zero: 



T <T r 



2(1 



So(q) 



p+q^o 



p+q|- 



(6.5) 



Now, taking into account the temperature dependence of the Bose-condensate fraction for 
an ideal gas [H H] 



no 
N 



T \3/2 



which is also determined from condition (5.10) for the temperatures T < T c . We find from 
(6.5) that 



p"|p+q| 



pt^o 
p+q^o 



In this expression we will pass from the summation over the wave vector p to the integration 
considering the thermodynamic limit iV — > oo, V — > oo, N/V — p — const and ultimately 
we find after integration over angles 



S (q) = coth I (3 



, k 2 q 2 
4m* 



2{T/T C 



\3/2 



+ 



m 



& _ i 4n 2 pqh 2 P 



P 



/3h 2 p 2 /2m* 



-In 



}H 2 (j>-q) 2 



dp, 



(6.6) 
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T < T r . 



For the temperatures which are above the critical temperature the Bose-condensate is 
absent (n = 0), that is why from (5.8) 



So(q) = 1 + 



to 



V 



x In 



An 2 pqh 2 P 



1 — Z e P 2m* 



-1 



e f3h 2 p 2 /2m* 



z e 



R h 2 (p-q) 2 
-P 2m* 



dp, 



(6.7) 



T>T C . 



We carry out analogous calculations also for the structure factor derivative (6.3): 



asb(g) 

dp 



T \ 3/2 



h 2 q 2 /Am** 



a^3L -i 

e P 2m* — 1 



TO 



lGir 2 pq(3m* 



f 

Jo 



dp ■ 



p° 



s 



3/2' 



sinh 



T < T r 



In 



1-e" 



1 - e _/ * 2m* 



(6.8) 



and 



dS (q) 
dp 



TO 



X 



1 - 



\Qn 2 pqPm 
d\nzo I h 2 p 2 



dp 



2m* 



In 



sinh 2 

1 - ^e"' 3 



p° 



a h 2 p 2 _ In go 

4m* 9 



4m* 

2m* 



z e 



(6.9) 



T > T r 



Let us now analyze the expression for the energy E from (6.1) in the low-temperature 
region. If T < T c then Zq — 1 and the contribution of the third term in (6.1) equals zero and 
the last two terms cancel each other, and this can be seen from (6.6) and (6.8) at P — > oo, 
S (q) -> 1, dS (q)/dp -> and also sinh ~ as well as cosh 2 [pE(q)/2] ~ 

e 0E{q) _ Consequently, we obtain Bogoliubov's formula for the energy as a mean value of the 
energy of the non-interacting elementary excitations: 



E = E Q + 



E(q) 



q^O 



e PE(q) _ I ■ 



T -> 0. 



In the quasi-classical limit H 
E 



from (6.1) we can easily find: 



2 2V 2^1+^ 



(6.10) 
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This expression corresponds to the random phases approximation also for the partition 
function (5.19). When arriving at it we took into account that at h — > the value of E 
in (6.1) gets canceled together with the second member in the square brackets of the third 
term since at high temperatures the high values of q prove to be important when a q — > 1. 
The fourth term in expression (6.1) at K — > also tends to zero and from the last fifth 
term we find a contribution taking into account that So(q) — > 1 and dSo(q)/d/3 = as can 
be seen from (6.7) and (6.9). Finally, the first term in the square brackets of the second 
term of the energy E is the energy of the ideal gas (with the temperature dependence of 
the effective mass m* being taken into account) and in the classical limit it equals 3NT/2. 

Hence, both at high and low temperatures our theory yields correct limiting cases. That 
is why we can expect good results also in the intermediate area of temperatures. 



7 STRUCTURE FACTOR AND THE POTENTIAL 
ENERGY 

By definition the structure factor is the mean quadratic fluctuation of particles density: 



S(q) = 





Pq! 


2 Rn 


(x\ 


x)dx 


I Rn{x\ 


x) 


dx 



(7.1) 



For calculating this expression we will use the designation of the partiton function (3.2) 
with integration and variables p q (5.3) and we will write (7.1) as a functional derivative: 

S(q) = -2 5 -^, (7.2) 

where the quantity v* is defined by formula (5.18). By means of expression (5.13) we will 
find from (7.2) that 



S ^ l + 5 (g)(a g tanh[/3E(g)/2]-tanh[^V/4m*])' ^ 

At the temperature of absolute zero (T = K) when v* = (a q — 1), we have the 
expression which is yielded by Bogoliubov's theory [THj l2~K]: 

S(q) = -. (7.4) 

The classical limit h — ► in (7.3) also yields a well-known result for the structure factor in 
the random phase approximation: 

s(9) = TTJW, <7 ' 5) 

We have one more illustration of the consistency of our theory. 
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Using formula (7.3) we find the mean value of the potential energy (2.5) calculated with 
the full Hamiltonian (1.1): 



Thus, 



J R(x\x)$dx 
J R(x\x) dx 



N(N- 1) 

w 



-u 



N x ^ 



q^O 



(7.6) 



So(q) 



1 + S (q) (a, tanh \J3E(q)/2] - tanh [(3h 2 q 2 /Am*]) 



In its turn, this expression makes it possible to calculate the mean kinetic energy if as a 
difference of the full energy (6.1) and (7.5): 



K = E-($) H . 



(7.7) 



8 EFFECTIVE MASS 



As we have already mentioned in [19] a different method was used to obtain a formula for 
the density matrix R^( product of the density matrix of the interacting particles 

with the effective mass m* multiplied by the factor P^(x'\x) which takes into account the 
interparticle correlations. Moreover, both for P^{x'\x) and for m* an expression was found 
only in the limit T = K. The effective mass is determined by the following equation [IT 



m 
m* 



r 



1 (q 

3iV ^-i a q (a q + 1) 



(8.1) 



q^O 



It is curious that this result coincides with the expression for the effective mass M* of 
the impurity atom moving in the liquid 4 He which we obtained in [21 using a completely 
different method 



M 
~M* 



1 - 



IN Z-w 



V, 



3N ^ Ot q (1 + Oi q ^) 
q^O " \ 1m) 



(1_M) 

V m J 



« 9 M 
y m 



(8.2) 



should the atom of the impurity be taken as equal to the mass m, i. e., to the mass of the 4 He 
atom. Expression (8.1) with the square bracket replaced by one was obtained independently 
in several studies [3711 ITTl l32l 133] . 

The fact that the "analytical extension "of formula (8.2) in the case when the impurity 
mass equals that of the liquid atom takes us to expression (8.1) testifies to a certain self- 
consistency of the theory [19J. The same coincidence of effective masses as shown in [19] 
cannot be expected at the temperatures not equaling zero. A dependence between the 
effective mass of the impurity M* on the temperature was studied in [3i| . 
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It is obvious that for high temperatures the role of interatomic interaction decreases 
and m* — > m when (3 — > 0. Let us consider this question in more detail. Let the impurity 
atom of the mass M and the coordinate R move in the liquid. The complete Hamiltonian 
of the system "an atom plus liquid"consists of the sum of the Hamiltonians of the liquid 
(1.1) plus the Hamiltonian of the impurity 

Mv 2 N ,. r , x 

ff^y + ^o + y^ ^e iqR Pq . (8.3) 

q^O 

Here the first term is the kinetic energy of the impurity atom whose velocity equals v, the 
two other terms being its potential energy of the pair interaction with the liquid atoms 
with the Fourier v q coefficient. As we consider a classical system its partition factor Z and 
aditional energy AE to the full energy of the liquid (6.10) can be easily found in the random 
phase approximation (see (5.19)) 

Z = Z N Z t , (8.4) 

where 

3 / 2 (AT 1 (RKi 



and Zn can be found from formula (5.19); 




Let us now introduce some effective Hamiltonian of the system "an atom plus liquid "in 
which the mass of the impurity atom is substituted by the effective mass M*, the interaction 
between the liquid particles v q is substituted by v q at q ^ and the interaction between the 
impurity atom and the liquid is absent. Constants z/ and Vq remain invariable at q = as 
we are interested by the renormalization via fluctuation mechanism, when 0. We will 
introduce the quantities M*, v q in the way so that the thermodynamic function of the initial 
system "an atom plus field "and the approximating model with the effective Hamiltonian 
should coincide in the respective approximations. For the suggested model with the effective 
Hamiltonian the partition function can be found easily 



3/2 



z = 'aMtttt^) ^ 

+ 



( M 

q^O q^O 



a N(N-l) R N_ 

2V V 



(8.7) 



Let us now compare expressions for the partition function of the initial system(8.4) and 
(8.5) and the model system (8.7): 



3 . M lv^ 
- In > 

2 M* 2 ^ 

q^0 



N \ N 
ln( 1 + /3-PJ -Pyi> q 
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— Y- 

2N ^ 1 



A" \ N 



(8.8) 



Taking into account some arbitrariness in the designation of M* and v q we will equate 
in (8.8) logarithmic terms. As a result of this we will find the effective potential 



and in the linear approximation the effective mass: 

M 



i-^>Jtt^tI • (8-9) 



Let us now put in (8.9) that M = m, v q = v q and find an expression for the effective 
mass m* in the classical limit H — > 0: 



m 1 ( Py v i 



- = 1 -^It*,'- (810) 

Thus, formulae (8.1) and (8.10) determine the effective mass in the limiting cases T = 
and H — > 0. One can rewrite these two formulae using the liquid structure factor S(q). At 
the absolute zero temperature from (8.1) and (7.4) one can find 

and from (8.10), (7.5) the high-temperature limit is obtained: 



m 

m* 3N 



i-^E^)- 1 ! 2 ' T ^°°- ( 8 - 12 ) 

q^O 

The above expressions might be "sewed together" numerically and the temperature depen- 
dence of the effective mass can be obtained. However, both of them are from the source 
other than the presented paper. Although, this is necessary at T = as the mass m* at 
this condition just drops out from all the expressions, as seen from (5.1). 

To reach full self-consistency of our theory at T ^ we need to find the expression for 
m* within the considered approach. One can demand, e. g., by means of special selection of 
m*, the expressions for the kinetic energy (7.7) and that from (K) = —mdF/dm (following 
from the know theorem about the derivation of the free energy F by some parameter) or 
(K) calculated from the density matrix R(x'\x) from (5.1) to equal. It is clear that in the 
exact theory all these expressions coincide, but as a consequence of the pair correlations 
approximation we use here they differ. To reach the self-consistency of the theory one 
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can utilize the arbitrariness of the parameter m* and choose it every time to obtain these 
formulae to be equal for the kinetic energy at different temperatures. 

A simpler approach to fix m* also exists. One can, e. g., choose m* to agree formula 
(7.3) for the structure factor with its experimental values at different temperatures. To 
make it in average for all the values of the wave vector q one can rewrite expression (7.3) 
as follows: 



1 1 

— r^r - a q tanh 

So{q) 



S(q) 



(3E{q) 



tanh ( (3 



, h 2 q 2 
4m* 



Let us make the integration over q in this equation, adding previously unity in both sides 
for the sake of convergence. After the integration over angles we find: 



q 



1 



1 



S(q) S (q) 



a q tanh 



(3E{q) 



l\\dq 



q 2 ( 1 — tanh 



P 



h 2 q 2 
4m* 



dq. 



(8.13) 



Integral in the right-hand side is easily calculated and we finally arrive at: 



-— - a q tanh 
S{q) 



PE(q) 



q 



So(q) 



dq + 



m 



f3h? 



dq 

3/2 



V^F(V2-l)C(3/2) 



(8.14) 



This equation solves the problem of the effective mass m* calculation via known S(q) for 
non-zero temperatures. Let us note that according to (7.4) this quantity equals to structure 
factor at T = 0, and S (q) also depends on m*. Equation (8.14) can be used to "sew together" 
formulae (8.11), (8.12) being correct both at low and high temperatures. 

Let us indicate one more way to determine m* via the condition that the pair distribution 
function of 4 He at zero interatomic distance should equal to zero due to hard core in 
the interaction potential. Only the numerical calculation can show how all the described 
approaches to the determination of m* concord. 

Therefore, if the liquid structure factor measured experimentally [H3 EHO EZl QE| is used 
as a source information, the theory of liquid 4 He proposed in this paper gives the connection 
between the observable quantities in the whole temperature interval, including the vicinity 
of A-transition. Such an approach was utilized by the author in |391 H2] for the calculation 
of different physical properties of 4 He near the absolute zero, T — > K. 

I express my gratitude to my colleagues from the Department for Theoretical Physics of 
the Lviv University, especially to Andrij Rovenchak, for the discussions about the results. 
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